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Abstract 

A new method to calculate level densities for non-interacting Fermions within the constant-spacing 
model with a finite number of states is developed. We show that asymptotically (for large numbers of 
particles or holes) the densities have Gaussian form. We improve on the Gaussian distribution by using 
analytical expressions for moments higher than the second. Comparison with numerical results shows that 
the resulting sixth-moment approximation is excellent except near the boundaries of the spectra and works 
globally for all particle/hole numbers and all excitation energies. 
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I. PURPOSE 



Our interest in the dependence of various nuclear level densities on energy and particle number 
has been triggered byrecent experimental developments in laser physics. The Extreme Light 
Infrastructure (ELI) [1] will open new possibilities for extremely high-intensity laser interac- 
tions with fundamental quantum systems different from the traditionally considered atoms and 
molecules J2]. At the Nuclear Physics Pillar of ELI under construction in Romania, efforts are 
under way to generate a multi-MeV zeptosecond pulsed laser beam [3]. For medium-weight and 
heavy target nuclei interacting with such a beam, photon coherence can cause multiple photon 
absorption. With energies of several MeV per photon, the ensuing nuclear excitation energies may 
well amount to several 100 MeV. Depending on the time scale on which the excitation takes place 
and on the specific nucleon-nucleon interaction rates, collective excitations may be induced, or a 
compound nucleus be formed |4|]. A theoretical treatment of the latter process along the lines of 
precompound reaction models requires the knowledge of the total level density, of the densities of 
p-particle /i-hole states, and of the density of accessible states for particle/hole numbers and/or 
excitation energies that go far beyond what has been considered until now. That applies not only 
to the target nucleus but also to all daughter nuclei populated by induced particle emission during 
the interaction time of the laser pulse. 

pi 

The standard approach to level densities goes back to the pioneering work of Bethe [151] who 
calculated the total level density as a function of excitation energy with the help of the Darwin- 
Fowler method. Basically the same method was used in many of the later works D6|4li|] dealing 
with the density of j9-particle /i-hole states and related quantities. A beautiful review is given in 
Ref. [12]. The Darwin-Fowler method yields analytical expressions involving contour integrals. 
Their evaluation, although straightforward, becomes increasingly involved with increasing num 



bers of particles and holes and/or increasing excitation energy. The same is true for Refs. [113 



14] 



that account for the exclusion principle by explicit counting. Moreover, without explicit numerical 
calculation it is not possible within these approaches to establish general properties of particle- 
hole densities like the overall dependence on excitation e nerg y and/or particle-hole number. More 
recent works use a static-path approximation (Refs. Ml 5L ll 6M and papers cited therein) or account, 

in addition, for the residual interaction in an approximate way (Ref. [17] and references therein). 
J L 

The method of Ref. H 1 811 avoids contour integrals and determines (again numerically) the level 
densities directly as coefficients of polynomials. The order of these rises rapidly, too, with en- 
ergy and particle/hole number. In none of these approaches does it seem possible to deal with the 
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enormously large values of the various densities attained for medium-weight and heavy nuclei at 
excitation energies of several 100 MeV in a practicable way. That is why we develop a different 
approach in the present work. 

In this Letter we present an analytical approximation to the global dependence of partial and 
total level densities that takes full account of the exclusion principle, that is valid for a finite num- 
ber of single-particle states, and that holds for all excitation energies and particle/hole numbers. 
We prove analytically that the level density for particles or holes is for a constant-spacing model 
asymptotically Gaussian. We improve on the Gaussian using analytical results for the low mo- 
ments of the distribution higher than the second. Comparison with numerical results shows that 
the resulting sixth-moment approximation is very precise except near the boundaries of the spec- 
trum (where numerical evaluation is easy). Particle-hole densities follow by convolution. The 
attained analytical form of the global dependence of level densities on excitation energy and parti- 
cle number extends our understanding of characteristic nuclear properties into uncharted territory. 
Moreover, we expect our results to be an indispensable tool in the calculation of laser-induced 
nuclear reactions mentioned in the first paragraph. 

We calculate the various densities in the framework of a constant-spacing model for spinless 
non-interacting Fermions. To justify our choice we consider by way of example the partial level 
density p p (E, J, %) for p particles and p holes, a function of excitation energy E, total spin J, and 
parity n, for a system of non-interacting Fermions in three dimensions. For other densities the 
reasoning is the same. The partial level density is given by 0JJJ] 

2J+ 1 



Pp (E,J,7r) = (l/2)p p (E) 



2V2^al p 



The factor 1/2 accounts for parity. The last two terms of the product give the spin dependence, 
with a 2p the spin-cutoff factor. With spin and parity being accounted for, p p (E) is defined as 
the level density of spinless non-interacting Fermions that carry no angular momentum. We note 
that in preequilibrium theories, the interactions between Fermions neglected here are taken into 
account as agents for equilibration. In our model, the non-interacting Fermions are distributed 
over a set of single-particle states. Each subshell with spin j of the three-dimensional shell model 
contributes (2j + l) states to the set. For large excitation energy or particle-hole numbers, we must 
take into account the exclusion principle exactly. It is equally important to account for the finite 
binding energy of particles and for the finite size of the energy interval available for holes. Both 



3 



strongly affect the various level densities at large excitation energies. We do so using a single- 
particle model with a finite number of states. Moreover, we calculate the various densities using 
a constant-spacing model for the single-particle states. It is clear from the shell model that the 
model is not realistic at the high excitation energies of interest. Taking into account the multiplicity 
(2j + 1) of the subshells, we note that the single-particle level density of the shell model strongly 
increases with energy. We return to this point at the end of Section PVT. 

H. APPROACH 

We consider / spinless Fermions in a single-particle model with constant level spacing d and 
with a finite number u of bound single-particle states. In the ground state all single-particle states 
from the lowest (energy d) up to a maximum level (energy F = fd with F for Fermi energy) are 
occupied. The remaining b = B/d levels (with B for binding energy) are empty. Here / and b are 
integers, and we have u = f + b. Excited states are described as p-particle /i-hole states, with p 
counting the number of particles in single-particle states with energy larger than F and not larger 
than B+ F, and correspondingly h counting the number of holes with energy less than F. For non- 
closed shell compound nuclei and/or nuclear reactions induced by composite particles, the number 
of hole states h may differ from p. We calculate various many-body level densities for non- 
interacting particles: pn(p, E) is the level density versus energy E for p particles confined to an 
energy interval of length B, pp(h, E) is the level density for h holes confined to an energy interval 
of length F, Pfb{p, h, E) is the particle-hole state density defined analogously, and pu(A, E) is 
the total level density for A particles distributed over an energy interval of length U = F+B. With 
e = E/dand e integer we define the dimensionless density e) = Pb(p, E) d and analogously 
for Uf(h, e), Wfb(p, h, e), and co u (A, e). All densities denoted by cu are integers. 

We describe the method of calculation for co b (p,e), assuming for simplicity of notation that 
b is odd and shifting the energy such that the ground state of the p-particle system has energy 
(l/2)p(p + 1). The maximum energy is bp — (l/2)p(p — 1), and the center of the spectrum is at 



The level density u>b(p, e) is defined as the number of ways in which p Fermions can be distributed 
over the b available single-particle states such that the total energy equals e, i.e., 



(2) 




(3) 
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The calculation of co b {p, e) poses a purely combinatorial problem. With = (1/2) (6 — 1) we 
define new summation variables k\ = ni — (1/2) (b + 1), I = 1, 2, . . . , p that range from —(3 to 
With e' = s — e (p) that gives 

W6(p,e')= 4i+fe 2 +...+fcp,e' • ( 4 ) 

-/3<fci<fc 2 <...<A;p</3 

We determine cj&(p, e') in terms of its low moments. Changing the signs of all summation variables 
in Eq. © one can easily show that ui b (p, e') = ui b (p, —e') is even in e', so that all odd moments 
vanish. For the 2mth moment with m = 0,l,2,...we have 



m b (p,2m) = J2(e') 2m u b (p,e') 

e' 

= E 

-/3<k 1 <k 2 <...<k p <l3 



2m 



(5) 



Following Ref. OD we adopt an occupation-number representation for Fermionic many-body 
states. We represent each set {fc;} of integers in Eq. © as a 6-dimensional vector {ui, u 2 , . . . , v b } 
with entries Vj that take values zero and one. The set {k{\ is represented by choosing Vj = 1 
in the p positions k\ and zero otherwise. The sum over all {k{\ is replaced by the sum over all 
6-dimensional vectors, i.e., over all choices of Uj subject to the constraint J2j u j = V- Thus, 



m b (p,2m) = 5 P,vi+u2+...+v b I 

Vl,U2,...,V b ^ j 

Q2m 

' J p,U 1 +U2 + ... + U b 



2 m 



da 



2m 



E 4 



ex P{ "^J I/ j} 



(6) 



(7=0 



We multiply Eq. © with exp{pa}, sum over p, and carry out the summations over the Vj. This 
gives the partition function 



Z b (a, 2m) 







2 m 



P 



da 2m 



[ (1 + exp{a + aj}) 



j=~P 



(7) 



(7=0 



The moment m b (p, 2m) is the coefficient multiplying exp{ap} in an expansion of Z b (a, 2m) in 
powers of exp{a}. For m = we find m b (p, 0) = the correct result. For m = 1,2 we obtain 



m b (p,2) 
m b (p,4) 



P 

E^ 2 

P 

E^ 4 

WE'* 

v i=-/3 



6-2 
p — 1 

6-4 
p — 1 

2 



6-4 
p- 2 



6-4 

p — 3 



6-4 
p-2 



(8) 



From Eqs. dH) we obtain the normalized moments 

M b (p,2 m ) = ^f^. (9) 
m b (p, 0) 

III. ASYMPTOTICALLY GAUSSIAN DISTRIBUTION 

Eqs. ([8]) suggest that asymptotically (6> 1,ji> 1) Ub(p, s') approaches a Gaussian distribu- 
tion. (Here, with 1 < p < b and b ^> 1, we consider p 3> 1 equivalent to p « 6/2. Particle-hole 
symmetry connects the cases p « 6 and p w 1). We recall that for a Gaussian distribution, the 
normalized fourth moment (see Eq. d9)) equals three times the square of the normalized second 
moment. For b 3> 1 and p 3> 1 that is exactly the relation implied by the values of m b (p, 4) and 
m b (p, 2) in Eq. ([8]). Indeed, taken by itself the last term in the expression for m b (p, 4) yields a 
value for M b (p, 4) which for 6 ^ l,p > 1 equals three times the square of M b (p, 2). More- 
over, the term proportional to j 4 in the expression for m b (p, 4) is smaller by the factor 1/p than 
the one proportional to (Yl J 2 ) 2 - To show that 0J b (p, e 1 ) becomes asymptotically (6 > l,p > 1) 
Gaussian we generalize the approach of Eqs. © to dH) to all moments. We define 

P 

G(a) = ' [ (1 + exp{« + aj}) = exp{H(a)} (10) 

j=-P 

and expand H(a) in a Taylor series around a = 0. With 

f(*) = j^h 01) 

1 + exp{a} 



and /w denoting the nth derivative of /, we have for m = 1,2,... 



<9cx r 

This shows that all odd derivatives of H vanish. We insert the Taylor expansion for H into Eq. ([101) 
and obtain 

r 00 i ^ i 

= G(0) exp £ a 2n f^ £ j 2 " . (13) 

^ n=l ^ j=-/3 ' 

From here, we proceed in two steps, (i) We neglect all terms with n > 1 on the right-hand side 
of Eq. (fT3l and show that as a result, co b is Gaussian for 6 ^> 1, p ^> 1. (ii) We show by complete 
induction that all terms with n > 1 in Eq. (fl~3l) are negligibly small in the same limit, (i) For 

co b (p, e') to be Gaussian we have to show that M b (p, 2m) = (2m — 1)!! [M b (p, 2)} m . Taking into 
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account the term with n — 1 only, expanding the exponential, and using the result in Eqs. (flOl) and 
© we obtain 

7=-j8 ' 

The normalized 2mth moment M b (p, 2m) is the coefficient multiplying expjjoa} in the expansion 
of Z b (a,2m) in powers of exp{a} divided by the normalization factor ( ). For b ^> 1 and 
p > 1 the relevant coefficient in G(0)(f') m ( b )~ is ~ \p(b - p)b~ 2 ] m . That yields M b (p, 2m) « 
(2m — 1)!! [M(,(p, 2)] m , consistent with a Gaussian form for co b (p,e'). In the last step of the 
argument we approximate products of the form p(p — 1) . . . (p — m) by p m . For fixed p the 
approximation becomes increasingly inaccurate as m increases. Our result is therefore valid only 
asymptotically, (ii) We use complete induction to show that the contributions of the terms with 
ri > 1 in Eq. (IT3T ) become vanishingly small for b 3> l,p ^> 1. We have shown above that 
this claim holds for n = 2 (i.e., for m b (p, 4)). We assume that the claim is correct for 2 < 
n < n , omit the corresponding terms in Eq. (Tl3T >. and show that it holds for n = n , i.e., for 
M b (p, 2n ). We have shown under (i) that the contribution to M b (p, 2n ) of the term with n — 1 is 
(2n - 1)!! (M b (p, 2)) no . The contribution of the term with n = n is G(0)f 2n °-^ £V j 2n °. For 
b > 1 we have £\ j m w 6 m+1 /(2 m (m + 1)). From Eq. (Qj) we have f = f - f. Therefore, 
y(2n -i) _ ^ 2 ™J C/ y« j s a polynomial of degree 2n in / with integer coefficients q, and the 
contribution of G{Q)f^ n °~^ J2j f n ° to M b (p, 2n ) is 

(by 1 b 2 ^ ^ /6-2n \ 
W (2n + l)2 2 ^^ Cl \ p-l J ■ 

For 6 ^> 2n and p ^> / we have ( 6 ~^™°) ~ {)p l /b l . The contribution (fT5l is, therefore, of 
order b 2no+1 while the contribution from the term with n = 1 is of order 6 3ri °. This shows that the 
contributions with n = 1 dominate all others. The situation differs for b ^> 1 and p w 1 or p 6 
where M b (p, 2) is of order 6 2 only and [M b (p, 2)} n ° is comparable in size to the contribution (fT5l ). 
Here the Gaussian approximation cannot be expected to work well. This is consistent with the fact 
that for p = 1 and p = b — 1 the densities are flat, co b (l,e') = 1 = u b (b — 1, e'). Furthermore, 
for p = 2 and p = b — 2 the densities have a triangularly shaped maximum. Only with p = 3 
and p = b — 3 does the density of states become Gaussian-shaped. The maximum at e' = (p) 
builds up only slowly as p increases from unity or decreases from 6—1. 
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IV. LOW-MOMENTS APPROXIMATION 



Using the asymptotically Gaussian form of cu b (p, e') we approximate that function in terms of 
its low even moments. We use Eqs. © for m b (p, 2) and m b (p, 4), calculate m b (p, 6) similarly, and 
find the parameters 7 2m , m = 1, 2, 3 of the normalized function 

ifW) = Cexp{- 72 (e') 2 - 74 (e') 4 -76(^) 6 } (16) 

that correspond to the normalized moments M&(p, 2m) with m = 1, 2, 3 in Eq. ©. The resulting 
function 



wr(p.O = y^ w (p,0 (17) 

is referred to in the following as the sixth-moment approximation to e'). Approximations 
obtained by using only the second (only the second and the fourth) moment(s) are denoted by 
ujf\p : e') (by o;^(p, e'), respectively). The same approach is used for Uf(h, e') and for u u (A, e). 
Except for suitable changes of indices and parameters, the results are formally identical. We 
mention in passing that the method is also useful for calculating the density of accessible states ifiol 



1911 under the constraints of the exclusion principle. Details will be given elsewhere. For the p- 
particle fo-hole density u>bf(p, h,e) we define e as the total excitation energy of the Fermionic 
system. Then 

u bf (p, h,e) = ^2 $e P +e h ,eUb(p, e p )u f {h, e h ) . (18) 



Here e p = e^°'(p) + e' p is the total energy of the p particles, and correspondingly for holes, while 
e = e p + e h is the excitation energy of the p-particle /i-hole system. Thus 

uj b fip, h,e) = Y^ ^ +£ / ft)£ _4°)_ £ £>)^&(P, £p)uf(h, e' h ) . (19) 

Since u b (p, sL) (u)f(h, e' h )) is a symmetric function of e' (of s' h , respectively), it follows that u b f 
is a symmetric function of e centered at = e p °^ + . Therefore we consider the function 
uibf(p,h,e') with e' = e — This function is symmetric about e' = 0. For the low even 
moments we obtain 



m bf (p, h,0) 

\pj \h / 

m bf (p,h,2) = m b (p,2) +m f (h,2) , 
m b f{p,h,A) = m&(p,4) + 2m b (p, 2)m f (h, 2) 

+m/(M), (20) 

and correspondingly for higher moments. 



V. NUMERICAL RESULTS 



We begin with an overview of the dependence of cub(p,e) on both p and e using the sixth- 
moment approximation (fTVT t . Even though we expect that approximation to work well only for 
p pa b/2, we display in Fig. Q] the values of usf\p, e) for b = 51 in the p-e plane as a coloured 
contour plot for all values of p between 3 and 6 — 3. For fixed p, the dimensionless energy e 
takes values in the interval (l/2)p(p + 1) < e < bp — (l/2)p(p — 1). This accounts for the two 
nearly parabolic and sawtooth-like boundaries of the coloured domain. The parabolic dependence 
on p is given by {l/2)p{p + 1) for the lower edge and by bp — (l/2)p 2 + {l/2)p for the upper 
edge. The contour plot is symmetric with respect to a simultaneous mirror reflection about the 
vertical line p = {b — l)/2 and about the horizontal line defined by the overall centroid energy 
e = (l/4)(6 2 — 1). This symmetry is due to the symmetry of Ub(p,s) in e about the centroid 
energy s^\p), and to particle-hole symmetry which equates u>b{p, e) with u>b(b — p, e) except for 
a shift by the difference (b — p)— (p) of the centroid energies. For fixed p, Ub(p, £) displays 
a maximum at ef \p) = {l/2)p{b + 1) (except for the cases p = 1 and p = b — 1 not displayed 
in the Figure). The location of the maximum increases linearly with p. This fact and the parabolic 
form of the boundaries cause the quasi-elliptical shape of the solid line of constant u^-values in 
the colour plot. We note the enormous maximum values of ujb{p, £ h (j>)) ~ 10 12 attained for 
p pa 25. All these features are generic (i.e., independent of the performance of the sixth-moment 
approximation) and apply likewise to uJf(h, s') and to u u (A, e'), except for a rescaling of abscissa, 
ordinate, and of the values of the densities. 

Limitations of the sixth-moment approximation become obvious when we consider the values 
of u)b(p, £) at the boundaries e = (\/2)p(p + 1) and e — bp — (l/2)p(p — 1) where we obviously 
must have Ub(p, s) pa 1. The sixth-moment approximation exceeds this value by one to two or- 
ders of magnitude, see Fig. |3] below. (We should keep in mind, of course, that the values at the 
boundaries predicted by the sixth-moment approximation are smaller by about 10 orders of mag- 
nitude than the values in the maximum. The relative accuracy of the sixth-moment approximation 
is, therefore, excellent). The white dashed lines in Fig. [TJshow at which values of p and e the 
sixth-moment approximation deviates by 10% from the exact values. The inaccuracy affects only 
the very tails of the density of states, symmetrically about the centroid energy e^ b °\p). The dashed 
lines thus follow the boundaries of uib(p, £)• 

We test the performance of the sixth-moment approximation in detail by a comparison with 
exact numerical results. This can be done throughout the critical domain (where either p pa 1 or 
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FIG. 1: Contour plot of the dimensionless level density oj^\p, e) for p particles in b = 51 equally spaced 
single-particle states based on the sixth-moment approximation (fTTT ) versus p and versus dimensionless 
energy e. The white dashed lines define the boundary of the region where w! 6 ' (p, e) deviates by 10% or 
more from the exact values. The full line presents the constant contour uyP'ijp, e) = 10 11 . 

p ~ b or where e is close to the boundary of the spectrum) since this is easily accessible numeri- 
cally. We calculate the exact values of e) in two ways, (i) We directly use Eq. ©. (ii) We 
use the occupation-number representation defined above Eq. © and sum over all 6-dimensional 
vectors {z/ 1; z/ 2 , . . . , z^}, grouping the results according to particle number p = £\ Vj and energy 
e' = Method (i) yields u)b(p,e) for fixed b and p. Method (ii) yields ujb(p,e') for fixed 

b and all values of p and e' . The demand on computing time is obviously larger for method (ii) . 



Particle-hole densities are then obtained from Eq. <TT81 > . Our results agree with those of Ref. 111011 
for the small numbers of particles and holes considered there. 

For the comparison between our exact and approximate results, we restrict ourselves to a few 
central features. In Fig. |2]we display the relative difference between the exact result and the sixth- 
moment approximation for b = 51 and various values of p near 6/2. Significant deviations occur 
only at the boundaries of the spectrum where the values of the sixth-moment approximation are too 
large. Even though the tails of the sixth-moment approximation are suppressed by many orders of 
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magnitude in comparison with the value at the center, that suppression is not strong enough. This 
is shown more clearly in Fig. |3] where we plot for b = 51 the values of oo^\p, e), of oo^ (p, e), 
and of wj® (p, e) (these functions are defined in and below Eq. (fTTT )) versus p at the boundary 
of the spectrum for b = 51. With every additional moment included in the approximation, the 
agreement with the correct value u;& ~ 1 is drastically improved. However, it would obviously 
take even higher moments than the sixth one to reach quantitative agreement at the boundary of the 
spectrum. This can be done, although convergence may be slow. Alternatively, we may calculate 
Ub(p, e) numerically for the critical values of p and e where the sixth-moment approximation is 
not sufficiently precise. These values lie at the boundaries of the spectrum shown in Figure [1] 
where either p fa 1 or p ^ 5 or e ^ (l/2)p(p + 1) or e ~ pb — {\/2)p{p — 1). In all these cases 
the number of terms that contribute to Eq. © is small, and the calculation is straightforward. 
Furthermore, for energies close to the spectrum boundaries , the density of states depends only on 
p and not on the number of levels b. For the lower boundary that is the case for [l/2)p{p + 1) < 
e < (l/2)p(p — 1) + b + 1. The numerical calculation can then be performed conveniently for 
a smaller number of particle states b' chosen such that Ub(p, e) = Ub>(p,e) for e in the energy 
interval of interest. For the dashed lines in Fig. [T] defining a 10% deviation of our approximate 
results, for instance, the exact values can be calculated numerically using b' ~ p(b — p)/10 + p— 1. 
It should also be borne in mind that in preequilibrium calculations one typically requires ratios 
(and not absolute values) of densities. We expect that these are predicted quite precisely by the 
sixth-moment approximation even at the boundaries of the spectrum. 

We turn to the total level density tu u (A,e) of A particles distributed over u equally spaced 
single-particle states as a function of the dimensionless excitation energy e. We have shown 
that u u (A,e) has approximately Gaussian shape, with a peak at half the total excitation energy 
(1/2) A(u + 1). The original calculation of cu u (A, e) by Bethe [5J] effectively also used a constant- 
spacing model but neglected the limitations due to finite particle number A and finite number u 
of single-particle states. With energy measured in units of the single-particle level spacing, the 
celebrated "Bethe formula" reads D5J 

WBethcOO = — \=exp{7i^2e/3} . (21) 
eV48 

We note that clothe (£) does not contain any adjustable parameters. The singularity at e = is 
due to the Darwin-Fowler method. The ensuing approximation fails at and near e = 0. Beyond 
this domain the Bethe formula yields a monotonically rising function of excitation energy since 
the underlying counting method assumes that the number of available single-particle states is 
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FIG. 2: Relative difference between the exact result and the sixth-moment approximation Alj^/u = 
[uj^\p, e') — Ub{p, s')]/uJb{p, e') versus energy e' for b = 51 and several values of p. The spectral bound- 
aries are indicated by the arrows on the abscissa. 

unbounded. Thus, there exists an excitation energy beyond which the Bethe formula exceeds 
co u (A,e) by an ever growing amount. This fact was qualitatively pointed out in Ref. 12011 . In 
Fig. ID we compare uiBetheiz) with the exact calculation for u = 51 and A = 41 and with the 
sixth-moment approximation for u = 250 and A = 200 (for the latter parameters the exact 
density of state values are not available for the entire energy spectrum). The latter parameter set 
mimics, very roughly, a heavy nucleus. Comparison with the exact calculation shows that WBcthc(^) 
underestimates the level density below the crossing point of both curves. We have found this to 
be a systematic trend. Both parts of Fig. |4] clearly display the crossing point and the increasing 
discrepancy between ^Beth e ( £ ) an d oj u (A, e) as e increases beyond this point. We interpret the data 
on the crossing points using the equilibrium distribution n(e) = (I /A) 1/(1 + exp{(e — A)/ (kT)} 
for A Fermions at temperature T with kT A and continuous single-particle energy e. With 
kT w y/e where e is the total excitation energy of the many-body system, we find that the crossing 
points occur at an excitation energy e where the fraction of particles in states with energies > u 
is of the order of a few percent. This is physically plausible. In a heavy nucleus this criterion 
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FIG. 3: Second-, fourth-, and sixth-moment approximation to cjj(p, e) at the lower spectrum boundary 

= (l/2)p(p + 1) versus p for b = 51. 

corresponds to excitation energies around 200 MeV. 

With increasing excitation energy, the constant-spacing model becomes increasingly unreal- 
istic. Indeed, the standard value 111 ill d ~ 13/ A MeV for the average spacing of single-particle 
levels near the Fermi energy in medium-weight and heavy nuclei strongly underestimates the 
single-particle level spacing in low-lying shells. We recall that every subshell with spin j con- 
tributes (2j + 1) states to p p (E) in Eq. ©. As a consequence, the number of states available for 
high-energy hole formation is smaller than predicted by the constant-spacing model. Therefore, 
the actual level density bends over more strongly with increasing excitation energy and terminates 
at a lower maximum energy than shown in Fig. |4] and the discrepancy with the Bethe formula 
is even bigger than presented there. The effect of an energy-dependent single-particle level den- 



sity was previously addressed, for instance, in Refs. [fell 12211 albeit under neglect of the exclusion 
principle. 

To account for the shortcoming of the constant-spacing model we are in the process of improv- 
ing our approach to calculate nuclear level densities. We divide the energy interval U = B + F 
into several sections / = 1,2,... with constant level spacing d t each but with d t ^ d v for I ^ I'. 
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FIG. 4: (a) Comparison of the exact level density oj u {A, e) for u = 51, A = 41 versus energy e with the 
Bethe formula (|21"I ). (b) The same for the sixth-moment approximation and u = 250, A = 200. 

Distributing p particles in all possible ways over these sections, so that there are pi particles in 
section /, we can use our results for the constant-spacing model in each section separately. For 
a fixed distribution {pi} the level density is a convolution over a product of Gaussians. The total 
level density is the sum over all distributions {pi}. We note that it will no longer be a symmetric 
function of energy. 
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VI. CONCLUSIONS 



Combining analytical and numerical methods, we have used a constant-spacing model for non- 
interacting spinless Fermions to develop a global approach to various nuclear level densities. This 
approach is viable also for large particle numbers and/or excitation energies where previous work 
runs into difficulties. As representative example we have displayed in detail the calculation of the 
particle density U(,(p, e) as a function of particle number p and excitation energy e. This function 
is symmetric about the center of the spectrum and, except for p = 1 and p = b — 1, displays a 
maximum at the center. It also possesses particle-hole symmetry. The shape of the boundaries 
of the spectrum and the two symmetries are responsible for the quasi-elliptical shape of the line 
of constant density in Fig. [TJ With Ub(p, e) ~ 1 at the boundaries of the spectrum, the value of 
tUb(p, s) at the center increases dramatically with increasing b and p « b/2, reaching values near 
10 12 already for b ~ 50 (and even larger values as b is further increased). The decrease by 12 or 
more orders of magnitude from the center of the spectrum to the boundary poses a considerable 
challenge to viable analytical approximations. Guided by the fact that for b 3> l,p 3> 1 cub(p, z) 
becomes asymptotically a Gaussian function of energy, we have used analytical expressions for 
the low moments to determine a sixth-moment approximation to e). This approximation is 
excellent except for values of p and e near the boundaries of the spectrum. These are indicated 
by the white dashed lines in Fig. [TJ Here the numerical calculation of Ub(p, e) based on Eq. © 
is easy and fast. Combining both approaches we obtain a reliable and easy-to-handle method of 
calculating the overall nuclear level density and p-particle /i-hole densities for medium-weight 
and heavy nuclei for all particle numbers and at all excitation energies. The results should be 
realistic except for limitations due to the underlying constant-spacing model. Because of shell 
effects the density of single-particle levels increases towards the Fermi energy, and this fact is 
not taken into account in the model. Work on a suitable generalization is under way. The Bethe 
formula is seen to fail beyond an excitation energy that amounts to approximately 200 MeV in 
heavy nuclei. 
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